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Abstract 

Contemporary studies in forensic entomology exhaustively evaluate gene sequences because 
these constitute the fastest and most accurate method of species identification. For this purpose 
single gene segments, cytochrome oxidase subunit I (COI) in particular, are commonly used. 
However, the limitation of such sequences in identification, especially of closely related species 
and populations, demand a multi-gene approach. But this raises the question of which group of 
genes can best fulfill the identification task? In this context the utility of five gene segments was 
explored among blowfly species from two distinct geographic regions, China and Pakistan. COI, 
cytochrome b (CYTB), NADH dehydrogenase 5 (ND5), nuclear internal transcribed spacers 
(ITS1 and ITS2), were sequenced for eight blowfly species including Chrysomya megacephala F. 
(Diptera: Calliphoidae), Ch. pinguis Walker, Lucilia sericata Meigen, L. porphyrina Walker, L. 
illustris Meigen, Hemipyrellia ligurriens Wiedemann, Aldrichina grahami Aldrich, and the 
housefly, Musca domestica L. (Muscidae), from Hangzhou, China; while COI, CYTB, and ITS2 
were sequenced for four species, i.e. Ch. megacephala, Ch. rufifacies, L. cuprina, and the flesh 
fly, Sarcophaga albiceps Meigen (Sarcophagidae), from Dera Ismail Khan Pakistan. The results 
demonstrate a universal utility of these gene segments in the molecular identification of flies of 
forensic importance. 
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Introduction 



Blowflies are the first insects to arrive at the 

scene of death, often seen ovipositing on the 

cadaver during the first few hours after death 

(Smith 1986). The growth periods, 

developmental rates (van Laerhoven 2008), 

and diapause responses (Ash and Greenberg 

1975) are all substantially different for closely 

related species. These peculiar features make 

blowflies the primary and most accurate 

indicators of the post mortem interval 

(Grassberger et al. 2003). However, the 

immature stages pose a great identification 

challenge for forensic entomologists due to 

their lack of species-specific anatomical 

characters (McDonagh et al. 2009). 

Identification of these stages using traditional 

methods (Mendonca et al. 2008) and/ or 

advanced techniques (McDonagh et al. 2009) 

is therefore the focus of many forensic 

studies. The most popular method of 

identification in recent years includes 

molecular taxonomy (Harvey et al. 2003; 

Smith and Baker 2008). In this regard several 

loci have been explored for their phylogenetic 

utility, but mitochondrial cytochrome oxidase 

subunit I (COI) (Harvey et al. 2003; Ratcliffe 

et al. 2003; Nelson et al. 2007; Ying et al. 

2007; Smith and Baker 2008) is often used, 

sometimes along with subunit II (COII) 

(Ratcliffe et al. 2003; Ying et al.,2007). The 5' 
■ . 

Table I. Species collected from Hangzhou, China and Dera Ismail 
Khan, Pakistan 



Species 


Family 


Location 


Chrysomya megacephaia (Fabricius) 


Calliphoridae 


Hangzhou, China 


Chrysomya pinguis (Walker) 


Calliphoridae 


Hangzhou, China 


Aldrichina grahami (Aldrich) 


Calliphoridae 


Hangzhou, China 


Lucilia sericata (Meigen) 


Calliphoridae 


Hangzhou, China 


Lucilia porphyrina (Walker) 


Calliphoridae 


Hangzhou, China 


Lucilia illustris (Meigen) 


Calliphoridae 


Hangzhou, China 


Hemipyrellia ligurriens (Wiedemann) 


Calliphoridae 


Hangzhou, China 


Musca domestica (Linnaeus) 


Muscidae 


Hangzhou, China 


Chrysomya megacephaia (Fabricius) 


Calliphoridae 


Dera Ismail Khan, Pakistan 


Chrysomya rufifacies (Macquart) 


Calliphoridae 


Dera Ismail Khan, Pakistan 


Lucilia cuprina (Wiedemann) 


Calliphoridae 


Dera Ismail Khan, Pakistan 


Sarcophaga albiceps Meigen 


Sarcophagidae 


Dera Ismail Khan, Pakistan 
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end of COI is also the site of the proposed 
universal animal DNA "barcode" (Hebert et 
al. 2003). COI barcodes have been 
successfully utilized in the identification of 
many blowfly species (Nelson et al. 2007). 
However, the barcoding approach has its own 
limitations (Moritz and Cicero 2004) as is 
seen with its failure to identify some closely 
related species of blowflies (Nelson et al. 
2007; Whitworth et al. 2007). Thus, COI 
barcodes, or for that matter any single gene, 
seems unlikely to resolve the identities of all 
calliphorid species of forensic importance 
because such phylogenies only infer the 
evolutionary relationships for the particular 
gene used. Consequently, a switch to multi- 
gene approach becomes necessary 
(McDonagh et al. 2009). To date only a 
handful studies have utilized such an approach 
(Wallman et al. 2005). 

Therefore, a multi-gene approach was 
employed for the identification of blowfly 
species of forensic importance from 
Hangzhou, China and Dera Ismail Khan, 
Pakistan. Gene data was freshly generated for 
a novel combination of five loci, including the 
barcode region of mitochondrial cytochrome 
oxidase subunit I (COI) (Hebert et al. 2003), 
NADH dehydrogenase 5 (ND5), cytochrome 
b (CYTB), and nuclear internal transcribed 
spacers 1 and 2 (ITS1 and ITS2). 

Materials and Methods 

Sample collection and Identification 

Samples of blowflies were collected from 
Huajiachi, Hangzhou over the period March 
2007 to June 2008, while in Dera Ismail Khan 
sampling was carried out in December 2008. 
Identification was performed using 
morphological keys (Fan 1992; Fan 1997; 



rnal of Insect Science.Vol. I I | Article 59 Zaidi et al. 

(■ \ 

Table 2. New sequences of flies from Hangzhou, China and Dera Ismail Khan, Pakistan 





GenBank accession numbers 


Species/individuals 


■TS, 


CO, 


CYTB 


ITS2 


ND5 


Hangzhou, China 


Chrysomya megacephalo 1 


FJ 7907 79 


CI / 1 A O 1 O 

FJ6 1 48 1 8 


FJ6 14835 


FJ6 1 4850 


CI t 1 Jfiffi 

FJ6 1 4869 


Chrysomya megacephalo! 


FJ 790780 


FJ6 1 48 1 7 


FJ6 14834 


FJ6 1 485 1 


FJ6 1 4868 


Chrysomya megacephalal 


rJ/VU/o 1 


C 1 i 1 A O 1 Z 

FJ6 ! 4o 1 6 


FJ6 1 4833 


ciz i ion 
FJ6 1 4852 


FJ6 1 4867 


Chrysomya pinguis 1 




FJ6 14820 


FJ6 14837 


FJ6 14853 


FJ6 1 487 1 


Chrysomya pinguis! 




FJ6I48I9 


FJ6 14836 


FJ6 14854 


FJ6 14870 


Hemipyrellia ligurriens 1 


FJ790782 


FJ6 14823 


FJ6 14840 


FJ6 14856 


FJ6 14874 


Hemipyrettia ligurriens! 




FJ6 14822 


FJ6 14839 


FJ6 14855 


FJ6 14873 


Hemipyrellia ligurriens! 




FJ6 14821 


FJ6 14838 


FJ6 14857 


FJ6 14872 


Lucilia sericata 1 




FJ6 14824 


FJ6 14841 


FJ6 14858 


FJ6 14875 


Lucilia sericata! 




FJ6 14825 


FJ6 14842 


FJ6 14859 


FJ6 14876 


Lualia sericata 3 




FJ6 14826 


Fj6 14843 


FJ6 14860 


FJ6 14877 


Lucilia illustris 1 


FJ790783 


Fj6 14827 


F]6 14844 


FJ6 14861 


FJ6 14878 


Lucilia illustris! 




FJ6 14828 


FJ6 14845 


FJ6 14862 


FJ6 14879 


Lucilia porphyrina 1 


FJ790784 


FJ6 14829 


FJ6 14846 


FJ6 14863 


FJ6 14880 


Lucilia porphyrina! 




FJ6 14830 


FJ6 14847 


FJ6 14864 


FJ6 14881 


Aldrichina grahami 1 


FJ790785 


FJ6 14831 


FJ6 14848 


FJ6 14866 


FJ6 14882 


Aldrichina grahami! 




FJ6 14832 


FJ6 14849 


FJ6 14865 


FJ6 14883 


Dera Ismail Khan, Pakistan 


Sarcophaga albiceps 




GQ9 12668 


GQ9 12671 


GQ9 12697 




Lucilia cuphna 




GQ9 12667 


GQ9 12670 


GQ9 12696 




Chrysomya megacephalo 




GQ9 12669 


GQ9 12672 


GQ9 12698 




Chrysomya rufifacies 1 




GQ9 12645 


GQ9 12673 


GQ9 12674 





Senior- White et al. 1940). The identified 
species are listed in Table 1 . 

DNA extraction, amplification and 
sequencing 

Twenty-six individuals including 22 blowflies 
(7 species) from Hangzhou and 4 individuals 
(including 4 species) from Dera Ismail Khan 
(Table 1) were used for DNA extraction, 
amplification, and sequencing. Total DNA 
was extracted from leg and thorax regions of 
each adult fly using the DNeasy Tissue Kit 
(QIAGEN, www.qiagen.com ) following the 
manufacture's protocol. Previously reported 
primers along with protocols of polymerase 
chain reaction (PCR) conditions were used to 
amplify the COI, ITS2 (Nelson et al. 2007), 
ND5 (Zehner et al, 2004), ITS1, and CYTB 
(Kengne et al. 2007). Thermal cycling was 
performed in a Mastercycler (Eppendorf, 
www.eppendorf.com ). Amplified products 
were purified and sequenced by Shanghai 
Invitrogen ( www . invitro gen.com ). The newly 



sequenced blowfly specimens from Hangzhou 
and Dera Ismail Khan with their related 
GenBank accession numbers for ITS1, COI, 
ITS2, ND5, and CYTB are listed in Table 2. 
The previously published sequences utilized 
during the present study are presented in 
Table 3. 

Sequence analysis 

All sequences were partial except for ITS2. 
Each COI, CYTB, and ND5 sequence 
represents part of the corresponding genes of 



I N 

Table 3. Previously published COI sequences used in our 
phylogenetic analysis 



Species 


GenBank accession numbers 


Chrysomya pinguis 


EU66 1 327, EU66 1 325, EU66 1 326 


Chrysomya rufifacies (Macquart) 




Chrysomya megacephalo 


AJ42604 1 . EU66 1 323, EU66 1 322. EU66 1 32 1 


Hemipyrellia ligurriens 


EU66I328 


Hemipyrellia pulchra {Wiedemann) 




Lucilia porphyrina 


EU66 1 329, EU66 1 330, EU66 1 33 1 


Lucilia illustris 




Lualia bazini (Seguy) 




Lucilia sericata 


AJ4222I2 


Lucilia cuprina (Wiedemann) 




Aldrichina grahami 


EU66I320 


Calliphora vicina (Robineau-Desvoidy) 
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heuristic search options except that 100 
replicates of random stepwise additions were 
used. Models of DNA substitution were 
estimated in Modeltest 3.7 (Posada and 
Crandall 1998). For ML we used GTR + G 
model for COI, Cytb and ND5, and F81+G 
model for ITSII nucleotide sequences. 
Bootstrap proportions were obtained after 
1000 replicates by using 10 replicates of 
random stepwise additions of taxa. 

Results and Discussion 



M. domestica (GenBank accession 
EU 154477) with COI pertaining to positions 
1385-2040 (655 bp), CYTB to 10722-11192 
(470 bp), and ND5 to 6340-6785 (445 bp) 
while ITS-II with its 422 bp length was 
completely sequenced. ITS1 sequences were 
not used in the phylogenetic analysis due to 
the difficulties in alignment as well as an 
inadequate data set. Nevertheless, ITS1 was 
sequenced for most of the species from 
Hangzhou with the legible length ranging 
between 400 to 700 bp. 



The software Dnadist (Felsenstein 1989) was 
used to compute distance matrices for 
nucleotide sequences of different gene 
segments and graphs were plotted using the 
same data. 

Chromas Pro 34-Version 1.33 (available 
online, 

www.technelysium.com.au/ChromasPro.html 
was used for manual editing of sequences that 
were subsequently aligned in CLUSTAL X 
using default parameter settings. 

Phylogenetic analyses were performed using 
maximum parsimony (MP) with PAUP* 
4.0M0 (Swofford 2002) and maximum 
likelihood (ML) with PhyML (Guindon et al. 
2005). The MP analyses were run with default 



Table 4. Intra-specific ITS I, COI and CYTB divergence among 
Chiysomya megacephala individuals 



ITSI Sequence Divergence 


Species/ individuals 


Ch. megacephalat 


Ch. megacephala! 


Ch. megacephala! 


Chrysomya megacephala 1 


0 






Chrysomya megacephala! 


1.20% 


0 




Chrysomya megacephala! 


0.90% 


0.30% 


0 


COI Sequence Divergence 


Species/ individuals 


Ch. Megacephala 1 


Ch. Megacephala! 


Ch. Megacephala 3 


Chrysomya megacephala 1 


0 






Chrysomya megacephala! 


0 


0 




Chrysomya megacephala! 


0.16% 


0 16", 


0 


CYTB Sequence Divergence 


Species/ individuals 


Ch. Megacephala 1 


Ch. Megacephala! 


Ch. Megacephala 3 


Chrysomya megacephala 1 


0 






Chrysomya megacephala! 


0 


0 




Chrysomya megacephala! 


0.24% 


0.24% 


0 



Sequence data were generated for 12 species, 
including 10 blowflies, from two geographic 
zones of Pakistan and China. 

Sequence variability for ITS 

The transcribed spacer ITS2 has been 
sequenced recently and successfully utilized 
in the molecular identification of blowflies 
(Nelson et al. 2007; Song et al. 2008). 
However, ITS 1 has not been sequenced before 
for blowflies. This region was found to be 
highly variable and was amplified with great 
difficulty. The species L. sericata, L. cuprina, 
Ch. rufifacies, and S. albiceps posed 
difficulties in amplification. Additionally, the 
baseline noise in chromatograms made the 
ITSI sequences unreliable for Ch. pinguis. In 
contrast, ITS2 was sequenced without any 
problems. 

Intra-specific sequence variation 

The ITS 1 intra-specific sequence divergence 
was computed only for Ch. megacephala. The 
variation was in a reasonable range, i.e. 0.29 - 
1.16 %, and was comparable to the 
mitochondrial segments of COI and CYTB 
genes (Table 4). However, previous studies 
usually showed no intra-specific variation in 
the Drosophilidae for ITSI region (Baffi and 
Ceron 2002). 
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The ITS2 region showed no intra-specific 
variation for most of the blowfly species 
except for L. sericata. This result seems to be 
consistent with an earlier study that suggested 
that ITS2 cannot be utilized in differentiation 
of geographical populations of some blowfly 
species (Song et al. 2008). 

ITS1 may be a good tool for identification at 
the population level. The most important 
reason is the high mutation rate that cause 
deletions or insertions. A prominent indel 
characteristic (a 184 bp long deletion) was 
detected in one of the Ch. megacephala 
sequences (FJ790781). The size and location 
suggested a loop deletion. Loop deletions, 
when present in coding regions, might have 
deleterious effects on the organism like any 
other mutation (Stefanovic et al. 2007), but 
their exact impact on non-coding regions such 
as ITS is unknown. In any case, use of indels 
(i.e. insertions and or deletions) as genetic 
markers has been recommended in 
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Figure I. Molecular variation at five loci for Lucilia illustris and other 
blowfly species. Note: Lucilia illustris is highly divergent from its sister 
species Lucilia porphyrina (—) at ITS I locus. High quality figures are available 
online. 



phylogenetic studies of natural populations 
(Vdli et al. 2008). Thus a high expectancy of 
insertion/deletions among ITS1 sequences 
might make them useful in identification of 
geographic populations of blowflies of 
forensic importance. 

Inter-specific sequence variation 

Inter-specific variation was explored for a 
wide array of paired species (blowfly-blowfly; 
blowfly- fleshfly, and blowfly-housefly). Both 
of the ITS regions were found highly variable 
as compared to mtDNA. The ITS1 segment 
was abnormally divergent among all gene 
segments ( Supplement Table SI ), ranging 
between 7.3-7.8 % (Figure 1). The results of 
the present study are consistent with previous 
studies of Dipteran families, such as 
Tephritidae, which also revealed little 
evidence of similarity between ITS1 
sequences of species, especially those 
belonging to different genera (Douglas and 
Haymer 2001). The most likely reason seems 
to be the rapidly evolving and mutating nature 
of non-coding regions, such as ITS (Haymer 
1994). On the other hand, ITS1 showed little 
variation among members of the same species 
as shown by Ch. megacephala individuals 
(Table 4). These results also showed no 
overlapping between intra- and inter-specific 
variations for ITS1. The ITS1 intra-specific 
variation as recorded for Ch. megacephala 
during the present study was 0.24 -1.16 % 
(Table 4), while inter-specific variation 
between different species pairs was in the 
range of 8- 60 % ( Supplement Table SI ). 
However, intra- and inter-specific overlapping 
seems plausible for ITS2 when the smallest 
inter-specific ITS2 value of divergence is 
recorded at 1.2 % for Ch. megacephala and 
Ch. pinguis ( Supplement Table SI ). Similar 
overlapping is previously observed among 
blowflies (Song et al. 2008), though without 
any taxonomic impact in that study, 
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Figure 2. Intra-specific ND5 sequence variation between Hangzhou and other regions. Note: Transitions are the most 
prominent type of mutations among ND5 sequences. High quality figures are available online. 



nevertheless overlapping can negatively 
influence the identification process. However, 
at the family level ITS2 possessed strong 
resolution power as it showed higher values of 
variation between blowfly-flesh fly pairs 
(14.7-21.1 %) and blowfly-housefly pairs (30- 
35.2 %), as compared to blowfly-blowfly 
pairs ( Supplement Table SI ). This result is 
also consistent with previous studies (Song et 
al. 2008). Another important finding of the 
present study was abnormally high ITS1 
divergence values for sister species L. 
porphyrina and L. illustris, i.e. 22 % (Figure 
1). At other gene loci, including ITS2, the 
divergence values between these two species 
were comparatively low, but adequate for 
differentiation (4-6 %, Figure 1). 

Sequence variability for mt DNA 

Mitochondrial gene segments showed less 
sequence variability as compared to the ITS 
regions although sequence divergence values 
were almost comparable between CYTB and 
ITS2 (Figure 1). Insertion and deletions were 
comparatively fewer in mtDNA genes. 

Intra-specific sequence variation 

These results showed a mean of 0.15% COI 
intra-specific variation for Ch. megacephala. 



The results are similar to those for Chrysomya 
species of blowflies (0.097 %) reported by 
Nelson et al. (2007). The mt gene CYTB 
showed variation of 0.24 % among both Ch. 
megacephala and L. sericata individuals. On 
the other hand, ND5 sequences of blowflies 
were highly conserved with no intra-specific 
variation. However, comparisons with ND5 
sequences from other regions of China that 
recently became publicly available (Table 3) 
showed significant variation among the 
transition/transversion ratios and the indel 
characteristics between species from 
Hangzhou and their counterparts in other 
Chinese cities (Figure 2). These results show 
the presence of distinct geographic 
populations of blowflies or perhaps cryptic 
species. Mutations such as insertions, 
deletions, transitions, and transversions 
(Figure 2) are the result of, or an indicator of, 
divergence of natural populations, making 
them an important tool in the population 
identification studies (Vali et al. 2008). 

Inter-specific sequence variation 

The inter-specific variation for mtDNA 
segments, including COI, CYTB, and ND5, 
was found suitable for blowfly species 
diagnosis. However, for some sister species 
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S albicepsD 



Figure 3. Maximum-parsimony phylogram based on COI (A), CYTB (B), 
and ITS2 (C) sequences of blowflies from Dera Ismail Khan, Pakistan and 
Hangzhou, China. High quality figures are available online. 



this variation was considerably low, e.g. 0.78 
% of COI (L. sericata vs L. cuprina) and 0.5 
% of CYTB (L. sericata vs L. cuprina). Even 
ITS2 that provided a least variation of 3.5 % 
for the L. sericata and L. cuprina pair showed 
low variation for another blowfly species pair, 
i.e. Ch. megacephala and Ch. pinguis (1.2 %). 
These results justify the need of a novel 
molecular marker for sister species diagnosis. 

Phylogenetic analyses 

In an attempt to compare blowfly species from 
two distinct climatic zones, single (Figures 
3A, B, C) and multi-gene trees (Figures 4A, 
B; Figure 5) based on COI, CYTB, ND5, and 
ITS2 sequences were constructed utilizing 
maximum parsimony and maximum 
likelihood methods. Both approaches 
identified the 11 fly species of forensic 



S albicepsD 



Ch rufifaciesDI 



1.00 



Ch pinguisHI 



Ch pinguisH2 
Ch megacephalaH2 
* «n - Ch megacephalaDI 
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I 00 f * grahamim 
A grahamiH2 



1.00 r L 
0.971 L 
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H ligurriensH3 
H ligurriensH2 
L cuprinaD 
L sericataH2 
L sericataH! 
sericalaH3 
r LillustrisHI 
I L illustrisH2 

L porphyrinaHl 



1 0(l 



L porphyrinaH2 



Figure 4. Maximum parsimony (A) and maximum likelihood (B) 
phylograms based on COI, CYTB, and ITS2 sequences of blowflies from 
Dera Ismail Khan, Pakistan and Hangzhou, China. High quality figures are 
available online. 



importance from China and Pakistan. This is 
the first time that the molecular data from the 
regions of Dera Ismail Khan, Pakistan and 
Hangzhou, China are presented and compared 
for blowflies. 

Single vs. multi-gene approach 

The significance of single genes in molecular 
forensics is undeniable (Nelson et al. 2007; 
Wells et al. 2007; Song et al. 2008). 
Nonetheless, additional genes become 
essential in challenging identifications 
(Nelson et al. 2007). A reasonable group of 
genes not only clarify doubts and provides 
validation, but also shed light on new 
evolutionary relationships. A unique instance 
is that of A. grahami. Grouped together with 
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Figure 5. Maximum-likelihood phylogram based on COI, CYTB, ND5, 
and ITS2 sequences of blowflies. High quality figures are available online. 



Lucilini tribe in the COI gene tree (Figure 3A) 
this species presented a deviation from 
traditional taxonomy because this species 
(Calliphorinae) is identified as sisters with 
Chrysomyiinae rather than Luciliinae (Rognes 
1997). Contrary to COI, the CYTB and ITS2 
trees (Figures 3b, c) followed this 
classification. These two distinct patterns of 
evolution were also observed previously for 
blowfly phylogenies based on 28rRNA 
sequences (Stevens and Walls 2001). The 



multi-gene trees validated both views. The 4- 
gene ML tree (Figure 5) agreed with the 
classical taxonomy while the 3 -gene MP and 
ML trees (Figures 4A, B) displayed the 
alternative scheme. 

However, each of the single gene trees, 
although they adequately identified blowflies, 
they also showed some limitations in terms of 
evolutionary relationships between species. 
For example, the ITS2 tree incorrectly put H. 
ligurriens within genus Lucilia (Figure 3C). 
Previously another single gene phylogeny 
based on the COI locus also grouped this 
species with L. cuprina (Wells et al. 2007). 
The multi-gene trees made amendments by 
supporting the monophyletic status of genus 
Hemipyrellia (Figures 4A, B). Chrysomya 
rufifacies provided another example of 
misplacement by single gene tree, CYTB in 
this case (Figure 3B). The 3-gene ML tree 
clarified this ambiguity by rightly placing this 
species along with other species of genus 
Chrysomya (Figure 4B). 

These results demonstrate a universal utility 
of the respective five gene segments. Though 
single genes exhibited easy identification, the 
wrong placing of some species illustrates the 
possibility of future misdiagnosis. Therefore a 
switch to a multi-gene tactic was made with 
future work aiming at the examination of 
more variable segments and species from 
distinct and multiple geographic regions. 
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Supplement Table SI. Pair-wise comparison of inter-specific variation among different gene segments 
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